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ABSTRACT 

Diffuse emission is often challenging since it is undetectable by most instruments, which are generally 
dedicated to point-source studies. The ^®A1 emission is a good illustration: the only available ^®A1 
map to date has been released, more than fifteen years ago, thanks to the COMPTEL instrument. 
However, at the present time, the SPI spectrometer aboard the INTEGRAL mission offers a unique 
opportunity to enrich this first result. In this paper, 2 x 10® s of data accumulated between 2003 and 
2013 are used to perform a dedicated analysis, aiming to deeply investigate the spatial morphology 
of the ^®A1 emission. The data are first compared with several sky maps based on observations at 
various wavelengths to model the A1 distribution throughout the Galaxy. For most of the distribution 
models, the inner Galaxy flux is compatible with a value of 3.3x 10“^ ph cm“^ s“^ while the preferred 
template maps correspond to young stellar components such as core-collapse supernovae, Wolf-Rayet 
and massive AGB stars. To get more details about this emission, an image reconstruction is performed 
using an algorithm based on the maximum-entropy method. In addition to the inner Galaxy emission, 
several excesses suggest that some sites of emission are linked to the spiral arms structure. Lastly, an 
estimation of the ®®Fe line flux, assuming a spatial distribution similar to ^®A1 line emission, results in 
a ®®Fe to ^®A1 ratio around 0.14, which agrees with the most recent studies and with the SN explosion 
model predictions. 

Subject headings: Galaxy: general - Galaxy: structure - gamma rays: general - gamma rays: ISM - 
nuclear reactions, nucleosynthesis, abundances - methods: data analysis 


1. INTRODUCTION 

The 1.809 MeV line emission associated with the 
^®A1 decay is one of the most intense gamma-ray 
lines observed in our Galaxy. It was first detected 
by the Ge spectrome ter on the HEAO-G spacecraft 
(jMahonev et al.l 119841 1. However, to date, only the 
GOMPTE^ Compton telescope aboard the Compton 
Gamma Ray Observatory has mapped the ^®A1 during 
its nine years survey. The emission has been found 
mainly distributed alo ng the Galactic plane and supports 


a massive-stars origin 1 

Diehl et al.lll995l lOberlack et al.l 

119971: IKnodlseder et al. 

Il999al Pliischke et al.ll2001ll. In 


addition, the early COMPTEL skymaps suggest a num¬ 
ber of marginally significant spots, some of them be¬ 
ing poten tially associated w ith the Galactic spiral arms 
structure (|Chen et al.lll996ll . However, most of these fea- 
tures remain compatible w ith statistical noise in the data 
(jKnbdlseder et al1ll999bll . 

These COMPTEL maps have been used as a basis to fix 
the spatial morphology of the ^®A1 line emission for sub¬ 
sequent works related to the spectral analyzes. Among 
them, detailed studies in the inner Galaxy and extended 
regions along the Galactic plane made with INTEGRAL 
SPI indicate that the intrinsic line width is less than 
1.3 keV and that line position shifts along the plane 
corresponding to the rotation of our Galaxy, confirm¬ 
ing at least a partial association of the ^®A1 emission 
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^ The COMPTEL telescope JSchoenfelder et al.l 1199311 was 
studying the band 1-30 MeV, over a wide field-of-view of about 
1 steradians. Around 1.8 MeV, the energy resolution is about 140 
keV and angular resolution '^4°. 


with the spiral arms (iDiehl et al.ll2005 iWang et al.ll2009l: 

iKretschmer et al.ll2013L and references therein). 

A related topic is the ®®Fe emission, more precisely the 
isotope lines at 1.173 and 1.333 MeV released when it 
decays into ®®Co and ®®Ni, for the final stage of massive 
star evolution. Because of its weakness, the ®°Fe radia¬ 
tion has been detecte d from our G alaxy wit h only two in¬ 
struments, RHESS I ( SmithI [20041 1 and SPI (jHarris et al.l 
120051: IWang et al.l 120071) . The ®°Fe to ^®A1 flux ratio is 
found to be about 0.15 by both instruments. 

In this paper, 10 years of INTEGRAL observation are 
used to examine the spatial morphology of the ^®A1 line, 
through direct sky-imaging and sky distribution model 
comparison. We will emphasize on the data analysis, 
especially the instrumental background modeling issue, 
a key point for both COMPTEL and SPI instruments. 
We have developed several methods, and try to derive 
reliable conclusions, independent of any specific (sky or 
background) model. 

In the following sections, we review the main charac¬ 
teristics of the instrument and data-set and discuss the 
basic principles of the developed methods in Section [2j 
We then present the results (Section (S]) on the global 
morphology of the ^®A1 line emission, but also on spe¬ 
cific regions suspected to harbor A1 progenitors, before 
discussing the outcomes in Section IH 


2. DATA AND ANALYSIS METHOD 
2.1. Instrument and observations 

The International Gamma-Ray Astrophysics Labo¬ 
ratory {INTEGRAL) observatory was launched from 
Baikonour, Kazakhstan, on 2002 October 17. 

The on-board SPI spectrometer is equipped with an 
imaging system sensitive both to point-sources and ex- 
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tended/diffuse emission. It consists of a coded mask as¬ 
sociated with a 19 Ge detector camera. This leads to 
a spatial resolution of ^2.6° over a field-of-view of 30° 
(jVedrenne et al.ll2003fl . The instrument’s in-flight per¬ 
formance is described in lBooues et al.l (120031 1. Due to its 
non-conventional coded mask imaging system, the imag¬ 
ing capability relies also on a specific observational strat - 
egy based on a dithering procedure (iJensen et al.ll200§l . 
where the direction of pointing of each exposure is shifted 
from the previous one by ^2.2°. A revolution lasts 3 
days, the time that the spacecraft performs a large eccen¬ 
tric orbit, but half of a day of data is unusable because of 
the radiation belts crossing. It contains generally about 
100 exposures lasting approximately 45 minutes each. 
The present analysis is based on public data recorded 
with the SPI instrument from revolution 44 (2003 Febru¬ 
ary 23) to revolution 1287 (2013 April 28). Around 1 
MeV, high-energy particles saturate the electronics and 
can generate false triggers. Nonetheless, it is possible 
to analyze the signal in this energy range thanks to an¬ 
other electronic chain (via Pulse Shape Discriminators or 
PSDs) not affected by the saturation problem. Th e pro¬ 
cedure is explained in iJonrdain fc Roau^ (|2009l l. We 
use the events which trigger only one detector (single¬ 
events). Note that the events which hit successively two 
or more detectors (multiple-events, representing ~25% 
of the photons in this energy band) are not used in this 
work. After excluding data contaminated by solar flares 
or the radiation belts, it results in about 77 000 expo¬ 
sures and 2x10® s observation livetime. 

We performed our analyzes of the ^®A1 line emission in 
the 1805-1813 keV band, to take into account the Ger¬ 
manium energy resolution (FWHM of 2.9 keV at 1764 
keV), including its degradation between two consecutive 
annealings (^ 5%). At these energies, the gain calibra¬ 
tion (performed orbit-wise) accuracy is better than ±0.01 
keV. 

2.2. Data modeling 

The signal recorded by the SPI camera on the 19 Ge 
detectors is composed of contributions from each source 
(point-like or extended) present in the field-of-view, con¬ 
volved by the instrument aperture, plus the background. 
For extended/diffuse sources, we assume that their spa¬ 
tial distributions are given by an analytical function or 
an emission map fSection l2.3l) whose intensities are to be 
determined. For a given energy band and for Ng sources 
located in the field of view, the data n^p obtained dur¬ 
ing an exposure (pointing) p for the detector d can be 
expressed by the relation; 

’^dp — ^ ^ dd,(ipjSj ± bdp ± Cdp ( 1 ) 

1=1 

where sj is the intensity of source j. Rdp,j is the response 
of the instrument to source j, bdp the background, and 
edp the statistical noise (both for exposure p and detec¬ 
tor d). We assumed that there are Np exposures and Nd 
detectors. 

At the energies considered in this paper (E > 1 MeV), 
point-source emissions are weak and stable in time within 
the measurement uncertainties. Concerning the ex¬ 
tended/diffuse sources, they are not expected to vary. 


Thus, for a given set of Np exposures and Nd detec¬ 
tors, the system of equations, as formulated above, re¬ 
quires Np X Nd equations to solve for Ng Np x Nd 
unknowns. Hence, it is mandatory to reduce the number 
of unknowns. We will see below that the observed prop¬ 
erties of the background allow us to strongly decrease 
the corresponding number of free parameters. Finally, 
to determine the sky model parameters, we adjust the 
data through a multi-component fitting algorithm, based 
on the maximum likelihood statistics. Expected counts 
are obtained by convolving a sky model with the instru¬ 
ment response and then adding the background model. 
The resulting distribution is compared to the recorded 
data with free normalizations of both components. We 
used Poisson’s statistics to evaluate the adequacy of the 
various sky models to the data. The core algorithm de- 
veloped to handle suc h a large system is described in 
iBouchet et al. I (j2ni3all . 

2.3. Modeling of the Al spatial distribution 

A way to estimate the ^®A1 emission spatial distribu¬ 
tion over the Galaxy (i.e. the corresponding Sj term in 
Equation 1) is to represent it with some templates. In 
this stud y, we have used maps lis ted in Table [l] sim¬ 
ilarly to iKnbdlseder et al.l (|1999al l for the COMPTEL 
data. These maps, although the list is not exhaustive, 
emphasize some large-scale structures of the sky observed 
at particular wavelengths and associated to specific emis¬ 
sion mechanisms that we may relate to the Al emission. 
Specific treatments have been applied to some of them: 
The A[3.5/r] and A[4.9/r] maps are the NIR 3.5 and 4.9/r 
maps corrected for reddening using NIR 1.25/r map and 
averaging emission at latitude |6| > 4 0° to estimate 
the ze ro level emission as explained in iKrivonos et al.l 
(iroi and references therein). Note that, with this 
procedure, the extra-galactic component has been re¬ 
moved, but the resulting maps are not expected to 
have an accuracy better than ^ 1 0%. For the CO 
(|Dame. Hartmann &: ThaddeusI [200111 and the EGRET 
(> 100 MeV) (provided by the NASA/Goddard Space 
Fli ght Center) maps, we apply the pre-treatment detailed 
in IKnbdlseder et al.l ()1999aL and references therein). 
Hence, the peak emission around the Galactic Center of 
the CO map has been removed, while point-sources from 
the second EGRET catalog have been subtracted from 
the EGRET map, together with an isotropic intensity of 
1.5x 10“® ph cm“^ s“^ sr“^ to take into account the 
cosmic diffuse background radiation. 

We have introduced a parameter to quantify the differ¬ 
ence between all these maps, the contrast, which we 
defined as the ratio of the flux contained in the re¬ 
gion |Z| < 150°, |6| < 15° to the total flux enclosed in 
|1| < 180°, |6| < 90°. This ratio varies from 0.4, for the 
25/r MIR map, to almost 1, for the A[4.9/x] map. Eor 
reference, the EGRET map, which traces the interstel¬ 
lar gas/cosmic-ray emission, has a contrast value of 0.7. 
The maps up to EGRET (left part of the x-axis of Eig- 
ures [6] and IB are said to have low-contrast and those 
above (right part) high-contrast. 

2.4. Background determination 

The instrumental background corresponds to a more 
or less isotropic component due to particles hitting the 
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telescope or created inside its structure. It is the main 
contributor to the flux recorded by the detector plane 
and represents a key issue since the signal-to-noise ratios 
considered in this study are below 1%. 

However, we have to note that the background term in 
Equation (1), formally consisting of M = Nd x Np values 
(one per detector and per pointing), can be rewritten as: 

bdp — dp X Ud X tdp (2) 

Here u is a vector of Nd elements, representing the ’’uni¬ 
formity map” of the detector plane (background pat¬ 
tern), tdp the effective observation time for detector d 
and pointing p. The evolution of the background in¬ 
tensity is traced with dp, a scalar normalization coeffi¬ 
cient per pointing. Assuming that u is determined else¬ 
where, the number of unknowns (free parameters) in the 
set of Np X Nd equations is reduced to Ng -\- Np (for a 
background intensity varying on the exposure timescale). 
However, the background intensity {dp) does not neces¬ 
sarily vary so rapidly and the number of related back¬ 
ground parameters could be still decreased, according to 
the actual background evolution. 

We detail below both aspects (detector pattern and 
timescale evolution) of the background determination. 

2 . 4 . 1 . Background intensity variations 

Figured] shows the mean count rate (per detector) evo¬ 
lution with time. In our standard analysis, we generally 
assume that the background varies with a fixed timescale. 
We have tested several timescales from one exposure (~ 
2-3 ks) to one revolution (^ 2-3 days) and compared the 
results from a statistical point of view. Figure [2| shows 
the evolution of several indicators as a function of the 
chosen timescale. The ” reduced chi-square” is defined as 
x\/dof, where xi stands for the likelihood equivalent 
chi-squar43 and dof for the degrees of freedom. 

The reduced chi-square (or similarly xi — dof) is rather 
stable for timescales less than 12 hours; above, its value 
increases very quickly. As we aim to have as few as 
possible free parameters, in order to get a better con¬ 
ditioned system of equations, we can consider that a 
timescale of ^ 9 hours is a reasonable trade-off. In com- 
p arison, for the 505-51 6 keV band, we have concluded 
in iBouchet et al.l (|20ir)ll that the best quality and robust 
results are obtained with a value of ^ 6 hours. Note that 
the behavior of the above quantities does not depend on 
the assumed synthetic map nor on the pattern determi¬ 
nation method (Section I2.4.2|) . 

However, to impose segments of fixed duration is conve¬ 
nient and easy, but not necessarily the best way to pro¬ 
ceed since the background intensity can vary with vari¬ 
ous timescales along the mission. We have therefore ap¬ 
plied a segmentation algorithm develop ed to determine 
the va riability timescale of the sources (|Bouchet et al. I 
l2013bll to the background signal. We first subtracted the 

^ To obta i n a equivalent statistics for Poisson-distributed 
data, IMighelll 111993 1 proposes the usage of the maximum likelihood 
ratio statistic xf = 2 X et — nt + Ui X In^nijei), M is the 

number of data poi nts and e,; i s the model of the data n^. This 
is the Cash statistic IICashlll9731 for Poisson-distributed data (low 
counts); C = —2 X ln ( L ) = 2 X Ci — UiX In ( ei ), where L is the 

likelihood function, within the constant term ni—nt xln { ni ). 


source contribution from the total counts. A rough ap¬ 
proximation of the sky signal is enough for this purpose 
(its contribution to the total counts is weak). Then, for 
each exposure, counts are summed over all detectors, and 
for each revolution, the corresponding time-series (up to 
one hundred of exposures) is segmented. The number 
of segments and their lengths are adjusted in order to 
obtain the minimal number of segments ensuring x^/dof 
< I. Figure [3| illustrates the case where the background 
intensity varies with different timescales along the revo¬ 
lution. 

The background intensity has been found stable (fit with 
one segment) during 477 over 1081 revolutions. Finally, 
the total number of segments has been significantly re¬ 
duced since only ^3000 segments are required to describe 
the whole dataset, instead of the ^12 000 segments used 
when fixing a hours timescale. To compare qual¬ 
itatively these results with those obtained with fixed 
timescales, we plot in Figure [5] the statistical quantities 
corresponding to the background segmentation method 
with isolated symbols (arbitrary abscissa since no fixed 
timescale). All of them point toward an improvement 
of the fit quality when the background variability is de¬ 
termined with a flexible timescale. Finally, these (quasi 
model-independent) time segments have been used to de¬ 
scribe the background evolution in our subsequent ana¬ 
lyzes. 


2 . 4 . 2 . ’’Uniformity map” 

The uniformity map or background pattern {u in Equa¬ 
tion 2) can be fixed by hand before the fitting procedure 
by using “empty-field” observations, thought to contain 
no source signal. The dedicated SPI “empty-field” ob¬ 
servations are rare, but the exposures whose pointing 
latitude direction satisfies |6| > 30° constitute a good 
approximation since they contain only weak contribu¬ 
tions from sources, at the energies considered here. They 
amount to 20% of the observations and have been used 
for building a set of background templates for different 
periods along the mission (about one per 6 months). 

To quantify the properties of the background pattern for 
each period, we define the vector u as 


u{d) = 


Y.P dp{d) 
Y.ptp{d) 


for exposure p satisfying |b| > 30° 


However, in the case of diffuse emission, the high- 
latitude exposure fields may contain some signal and 
the background pattern deduced from them may be 
’’blurred”. More precisely, the high-latitude |6| > 30° re¬ 
gions contain ^ 30% to 40% of the total emission for the 
’’low-contrast” maps (12^, 25/r and synchrotron) while 
this ratio is less than 10% for ” high-contrast” (EGRET 
to A[4.9p,] maps) ones. Gonsequently, if the true emission 
distribution approaches the 25/i map, the detector pat¬ 
tern wifi be more affected by the high-latitude emission 
than if the true emission distribution follows the A[4.9 /j,] 
map. Note that the effect on the measured source flux 
is not predictable since the fit procedure adjusts back¬ 
ground and source normalizations: a “bad” background 
pattern may imply an underestimation or an overestima¬ 
tion of the source flux. 

In addition, for the ^®A1 study, we have to keep in mind 
that the side shields do not stop 100% of photons. This 
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means that any uniformity map or background pattern 
will contain also some diffuse emission signal passing 
through the shield. 

We have investigated another approach to estimate the 
detector uniformity pattern, by fitting it during the con¬ 
vergence procedure. We have to note, in this case, that a 
prerequisite is to have a sufficient knowledge of the source 
contribution and that the background model relies on 
more free parameters. With this alternative method, the 
Xl decreases by ^8200 for Adof=411 additional 
parameters (2973 parameters are used to describe the 
background). However, the improvement of the cri¬ 
terion is not reliable enough since the recovered source 
signal becomes background dependent and can be altered 
to an extent which is difficult to estimate. 

The two pattern determination methods are subse¬ 
quently referred to as fixed-pattern and fitted-pattern 
ones, and we have systematically compared the results 
obtained with each of them. 

Figure |4] presents the ^®A1 line flux obtained in the 
inner Galaxy (|1| < 30°, |6| < 10°) for both pattern de¬ 
termination methods and several background variability 
timescales with the same model for the ^®A1 distribution 
(here the 60^ template). With the fixed-pattern method, 
the recovered flux depends on the assumed background 
timescale, especially for background timescales below 
6 hours, while its value remains unaffected when the 
background pattern is adjusted. It is worth noting that 
the same analysis has been applied to the annihilation 
radiation, at 511 keV. Due to a better statistic and 
possibly a relatively less emission at high-latitudes, the 
fluxes obtained with the two pattern determination 
methods are perfectly compatible for each of the sky 
components, whatever the assumed background varia¬ 
tion time-scales. 


2.5. Imaging the sky 

Although the model-fitting process provides the best 
quantitative information about the global emission and is 
better suited to determine its level of confidence, a direct 
“imaging” algorithm provides more qualitative informa¬ 
tion such as the position of potential emitting sources 
and their extent as well as some basic, but model inde¬ 
pendent, characteristics of the emission morphology. 

To build an image, the sky is divided into small areas or 
pixels, the flux / in each pixel is to be determined. The 
linear model of the data derived from Equation (1) is put 
in a more synthetic form through 

n = Rf + b + e (4) 

where n represents the data, b the background, R the re¬ 
sponse of the instrument as described in Section [2. 2 1 and 
e the statistical noise, n, b and e are vectors of length M 
(the number of data points), / a vector of length N (the 
number of pixels in the sky) and R a matrix of size M by 
N. The statistical noise is assumed to follow a Gaussian 
distribution with a variance tr^ and null mean. 

However, when the number of pixels is large, the sys¬ 
tem of equations becomes ill-conditioned; the direct 
least-square solution is not always reliable and generally 
poorly informative. To select an informative and par¬ 
ticular solution, the most common technique consists in 


including a regularization term in the above system, in 
addition to the least square constraint. This solution, 
although biased by the choice of the regularization crite¬ 
rion, improves the conditioning of the system. A compro¬ 
mise between the goodness of the fit, quantified by 
and the regularization, say H, is found by maximizing 
the function 

Q(/)=ai7(/)-ix"(/) (5) 

where a is a parameter, which determines the degree of 
smoothing of the solution. 

We choose to use one of the most popular regularization 
operators: the entropy function. This metho d is known 
as the Maximum-Entropy Method or MEM (iGuHlllQ^ 
and references therein). Eor an assumed positive additive 
distributions, 

^ f- 

fdog^ ( 6 ) 

where mi is the default initial value assigned to the pixel 
i. Note that fi and mi are positive quantities. MEM has 
been proved to be a succ essful techniqu e for astronomi¬ 
cal image reconstruction (jSkillinglll981fl . The algorithm 
is described in Apppendix However, despite its ca¬ 
pabilities, it suffers from several shortcomings: among 
them, the difficulty of selecting the appropriate entropy 
function distribution and the default pixel values. In ad¬ 
dition, the possible correlations between pixels are not 
taken into account properly. One way to remedy this 
problem is described in Section 12.5.21 

2.5.1. Including the background intensity determination 

In our application, the instrumental background could 
be fi xed through the modeling method presented in Sec¬ 
tion However, such a modeling, even though be¬ 

ing accurate enough for the model-fitting procedure, 
may be a source of biases for the image reconstruc¬ 
tion, by preventing the appearance or by exaggerating 
the significance of some structures potentially present 
in the image. Consequently, we prefer to have the 
possibility to refine the parameters linked to the back¬ 
ground intensity during the image reconstruction pro¬ 
cess. The instrumental background variation timescale 
is fixed through the modeling method presented in Sec- 
tion l2.4l We search for both the solution vector / (length 
N) and the background intensities, b (length Nb). We as¬ 
sume that the background parameters are precisely de¬ 
termined/constrained (small error bars) compared to sky 
pixels, hence they do not require any smoothing. Note 
that the ” cross-talk” effects between sky pixels and back¬ 
ground parameters are reduced since the number of pa¬ 
rameters to describe the background is nearly minimum 
fSection l2.4.1ll . 

2.5.2. Pixel correlation 

As mentioned above, the possible correlations between 
the pixels of the image (the sky emission should vary 
smoothly from one pixel to the next ones) are not taken 
into account properly with the entropy function. It is 
possible to introduce artificially such a correl ation by us - 
ing an Intrinsic Gorrelation Eunction, ICE, (iGulllflOS^ . 
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The function to maximize becomes: 

Q(/) = aH{f) -\x\fiCF) where fwF = ICF{f) (7) 

For our work, we choose a radial Gaussian function. Such 
an IGF takes part largely to the solution regularization 
by ensuring its smoothness. 

2.5.3. Construction of a sky image 

In practice, the sky is divided into pixels of equal 
area, following the Healpix (Hierarchical E qual Area iso- 
Latit ude Pixelization of a sphere) scheme (iGorski et al.l 
1200511 . and initialized with a uniform default value. In 
a first step, we consider a very low-resolution map (48 
pixels of ~ 29° resolution) and select the value of a such 
that the reconstructed flux in the inner Galaxy is com¬ 
parable to the value obtained with the model-fitting pro¬ 
cedure. The size of this problem is relatively small and 
a N-dimensional search-algorithm is used to maximize 
the function Q (^3000 parameters). The calculation of 
the first-order solution (image and background param¬ 
eters) is based on a classical Newton type optimization 
algorithm with a positivity constraint on the solution. 
During this first stage, the background pattern is fixed. 
The resulting image is used as a template to compute an 
improved pattern which is then fixed for the rest of the 
process. We then use this solution as a starting point 
(for an initial guess of the background parameters, the 
image is assumed to be flat) to build a high-resolution 
image (49 152 pixels of 0.9° resolutio n) through an al- 

S "’m similar to that proposed by iSkilling &: BrvanI 
. The parameter a and the uniform sky default 
model are computed at each iteration so that the solu¬ 
tion is ensured to follow the optimum MEM trajectory 
(See Appendix lA|l . At the end, the final solution satisfies 
/dof ~ 1 and fulfills an additional test on the degree 
of non-parallelism between the gradient of C and H (see 
Appendix |A|. Note that, the greater the number of it¬ 
erations is, the lower the chi-square value but the image 
is more “spiky”. The above mentioned stopping criteria 
results probably in a ’’spiky” , but also more objective 
image (| Skilling fc Brvanlll984ll . 

3. RESULTS 
3.1. Imaging 

The image displayed in Figure [S] indicates that the 
emission is essentially confined in the inner Galaxy re¬ 
gion (|/| < 30°, |6| < 10°), with a flux in the inner Galaxy 
of 3.5xl0“'*ph cm“^ s“^ . Note that the morphology 
is drastically different from the electron-positron anni¬ 
hilation line emission one, which is essentially concen¬ 
trated in the bulge region (iWeidensnointner et al.ll2008t 
iBoiichet et al.l[20inl : iGhurazov et al.ll20ir)( l. 

Beyond the morphology of the diffuse emission, our inter¬ 
est has been piqued by several spots visible in the image. 
Simulations of image reconstruction (Appendix [B]) show 
that structures with a peak intensity at ^one-tenth of the 
image maximum intensity are probably due to residual 
statistical noise. Above this level, structures are worth 
being considered carefully. Indeed, some of them have 
positions compatible with sources detecte d or studied 
during recent SPI investigations (See section [3^ . For in¬ 
stance, the image reveals several excesses in the Gygnus 


region (72° < I < 96°, -7° < 6 < 7°). A significant one 
is located at a position co mpatible with Cyg O B2 clus¬ 
ter (81°,-1°) as reported bv iMartin et al.l (1200911 . In the 
Sco-Cen region (328° < I < 355°, 8° < b < 30°), the 
strongest excess is localized at (1, b) ~(360°, 16°), hav¬ 
ing a spatial extent with a radius of 5°. Some structures 
appear also around the Garina and Vela regions and two 
others, more extended, in the Taurus/Anticenter region 
(105° < I < 170°, -15° < b < 20°). Finally, an additional 
spot detected above Za is worth investigating. However, 
a model-fitting analysis is required to extract more quan¬ 
titative information on the above mentioned structures 
(Section 13.311 . 

3.2. Testing template maps 

This analys is is similar to that done with the COMP- 
TEL data by iKnodlseder et al.l (|1999all , and consists of 
using template maps to model the spatial distribution 
of the ^®A1 emission through the Galaxy. In this case, 
only the intensity normalization is adjusted (1 parame¬ 
ter) in addition to the background parameters. In prac¬ 
tice, we fit the data with a combination of one of the 
maps listed in Table [T] plus a background mo del. The 
template preparation is detailed in Section 12.31 
To qualify the detection significance of the ^®A1 emis¬ 
sion for a given template, we adjust two models to the 
data: the first one contains only the background while 
the second contains the background plus the tested tem¬ 
plate map. Quantitatively, the improvement of the like¬ 
lihood is 512 for 1 additional parameter for the 25p, map 
(best model) and the background determined with the 
fixed-pattern method. The ^®A1 flux detection signifi¬ 
cance is ^ 23(7, and the associated reduced-chi-square is 
Xh/dof = 1.0038. Indeed, each template leads to a sig¬ 
nificant flux detection. The worse case, H/ map, has a 
significance of ^ 17a. 

The same analysis with the fitted-pattern background 
determination gives an improvement of the likelihood of 
312 (for the 100/i map). The ^®A1 flux detection signifi¬ 
cance is ^ 18(7 and x\/dof — 0.9973. 

Figure [6] displays the Maximum Likelihood Ratio (MLR) 
obtained with both pattern determination methods 
(fixed and fitted-pattern) for each of the tested maps (or¬ 
dered by increasing contrast). Several maps give similar 
results, leading to the conclusion that star-related distri¬ 
butions (FIR and MIR maps, dust and free-free distribu¬ 
tions, all with rather high-contrast) give a good descrip¬ 
tion of the data, as expected from what we know about 
the ^®A1 emission process. Similarly, 25/i and 12/r maps 
constitute good tracers, while presenting a low-contrast. 
In another hand, H/ and 53 GHz synchrotron maps can 
be excluded. In fact, we have seen (Figure [S]) that the 
emission is confined into the central part of the Galaxy 
and close to the disk. This explains why the H/ map does 
not provide a good fit to the data since it extends far in 
longitude and contains important emission at |(| > 30°. 
We can also point out that distributions built by imaging 
method (MEM or MREM) from COMPTEL data appear 
as good tracers of the ^®A1 emission, as expected. 

As shown in Figure [71 the reconstructed flux attributed 
to the inner Galaxy does not depend much on the tem¬ 
plate map and is compatible with a value of 3.3 x 10“'* 
ph cm“^ s“* , in agreement with the previously re- 
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ported values. We observe that the flux is systematically 
higher by (low-contrast maps) to 20% (high-contrast 
maps) when the fixed-pattern background determination 
method is used, but it has no scientific implication. How¬ 
ever, the total flux integrated over the whole sky clearly 
depends on the contrast of the assumed model. The more 
the map is contrasted (from left to right on Figure [7]), 
the weaker is the reconstructed total flux. This is due to 
the fact that the recovered global intensity relies on the 
central parts of the image (both higher flux and signal 
to noise ratios) and that a contrasted map encompasses 
less flux in its external parts than a flatter map. More¬ 
over, we are aware that low-contrast maps may suffer 
from significant “cross-talk” between the low-spatial fre¬ 
quency structure (a kind of pedestal which mimics a flat 
low-surface brightness emission) and the (more or less 
uniform) background contribution. 

As an additional test, we have performed correlations be¬ 
tween the direct image with 6° resolution (e.g. ICF of 6° 
FWHM) and each of the templates downgraded to the 
same ^6° spatial resolution (except the DMR/COBE 53 
Ghz one, which originally has 7° resolution). Indeed, the 
linear correlation coefficient does not depend much on 
the template map and keeps a value greater than 0.9 in 
latitude, and above 0.7 in longitude, except the H/ map 
(coefficient of 0.4 in longitude). 

Finally, it appears that it is hard to firmly conclude about 
a unique solution. This reflects the difficulty of determin¬ 
ing precisely such an extended weakly emitting structure 
and the similarity presented by most of the considered 
maps. However, we note that the curves follow the 
same evolution regardless the background determination 
method and hence that the conclusions do not depend 
on it. 

3.3. Regions of potential excesses 

To check quantitatively the significance of the most 
significant excesses and known ^®A1 emitting regions, we 
have performed a more complete model-fitting analysis. 
Note, however, that our analysis is not optimized for ex¬ 
tended point-sources since a map contribution is neces¬ 
sarily subtracted at the position of the sources and, at 
worse, may make them disappear. 

The sky model consists of one of the templates listed 
in Table [H to which is added a spatial model includ¬ 
ing the spots. For the sources already detected or 
investigated at 1.8 MeV, we have used the positions 
and spatial extens ions provided by previous works (Vel a 
(jDiehl et al.||199 5ll. Cvgnu s region (jMartin et al.l 1200911. 
Sco-C en (jDiehl et al.ll201flll. Orion/ Eridanus. (iVoss et al.l 
l20ir)H and Carina (jVoss et al.ll201^ . indicated in bold in 
Table [T|). To investigate the additional spots not yet 
referenced as ^®A1 emitters, the extent and location are 
based on the image analysis using simple models (point- 
source, Gaussian or disk). However, given their low- 
significance and the SPI spatial resolution of 2.6 °, their 
identification with known sources is just indicative. We 
report in Table 2, the flux values obtained using the IC 
(low-contrast) and A[4.9/x] (high-contrast) maps, as rep¬ 
resentative of the global results. 

The Cyg OB2 cluster (81°,-1°) has a flux of (4.1 ± 1.5) x 
10“® ph cm“^ s“^ similar to the value of (3 .9 ± 1.1) x 
10“® ph cm“^ s“^ obtained bv iMartin et alJ (|2009ll and 


the (3 .7±l.l)x 10~^ ph cm ~^ s“^ obtained with COMP- 
TEL (iPliischke et al.l^OOOll . We note also another excess 

at (l,b)~(100°,6°) with a flux of (2.7 ± 1.1)x 10“® ph 

— 2 —1 
cm s . 

Str ong emission has b een reported in the Sco-Cen region 
by iDiehl et Ml (j201fll l centered around (1, b) = (350°, 
20°). At this position, we find a flux (4.1 ± 1.6)x 10“® 
ph cm“^ s“^ when using an Hj map to model the large 
scale ^®A1 emission. However in our image, the local 
maximum in this region appears shifted to (1, b) ~(360°, 
16°), and the source is less e xtended (disk radi us of 5° 
instead of 10° as reported bv iDiehl et al.l (|2010D b while 
the flux remains at (4.1 ± 1.1) x 10“® ph cm“^ s“^ using 
again the Hj map. On the another hand, in both cases 
the flux clearly depends on the template map used and 
is less than 2 (t for most of the cases. 

The Carina and Vela regions are marginally detected 
{2a). For Carina, the measured flux of (2.8 ± 1.5)x 
10“® ph cm“^ s“^ is comparable to the est imation of 
(1.5 ± 1.0) X 10“® ph cm“^ s“^ reported bv IVoss et al.l 
(|2012ll using SPI data, and (3.1 = 1:0.8) x 10“® ph cm“^ s“^ 
reported by iKnodlseder et al.l (1199611 using COMPTEL 
data. For Vela, the flux of (3.3 ± 1.8) x 10“^ ph cm“^ 
s“^ is comparable to the value (3.6±1.2)x 10“® ph cm“^ 
s“^ reported by IDiehl et ^ (1199511 using COMPTEL 
data. For the Orion-Eridanus area (180° < I < 210°, 
-30° < 6 < 5°), we get an upper limit for the total flux of 
3x 10“® ph cm“^ s“^ (1 a) in agreement with the value 
of (4.5 ± 2.1) X 10“^ ph cm“^ s“^ obtained bv iVoss et al.l 
(1201011 from a synthesis model based on the Orion star 
populations. We also mention the the 2g upper limit 
of 1.7 x 10“® ph cm“^ s“^ reported bv lOberlack et all 
(jl995ll for the Orion region using COMPTEL instrument. 
Concerning the spots not yet reported, we first point out 
two extended structures in the Taurus/Anticenter region 
(105° < I < 170°, -15° <b< 20°) with fluxes of (5 ± 2) 
and (8 ± 3)x 10“® ph cm“^ s“^. Finally, a spot worth 
mentioning since it is detected above 3 cr, it is located 
at high-latitude (l,b)Ri(226°76°) with a flux of (7 ± 2)x 
10“® ph cm“^ s“^. We did not find any convincing po¬ 
tential counterpart, but consider the stability of the at¬ 
tributed flux against the assumed ^®A1 diffuse emission 
map as a good criterion to assess that an excess is robust 
and reliable. 

We did not notice any systematic effect in the source flux 
determination between the fixed and fitted background 
pattern adjustment method. Note that adding all these 
excesses into the model of the sky improves the likelihood 
only marginally (Ay|^ ^ 5) compared to the case where 
they are neglected, and that dedicated analyzes must be 
conducted to refine individually each result. 

3.4. and to ratio 

While the intensity of the ®°Fe isotope emissions at 
1.173 and 1.333 MeV contains important complementary 
information for the star evolution study, the weakness 
of the flux makes it impossible to derive any constraint 
on the spatial distribution. The latter is, therefore, as¬ 
sumed to be the same as the ^®A1 line one lSection l3.2|) . 
which is physically reasonable since Al and ®°Fe are be¬ 
lieved to be produced at least part ly in the same s ites (e . 
g., massive stars and supernovae. iTimmes et al.l (1199511 . 

















































Galactic ^®A1 emission with INTEGRAL SPI 


7 


iLimongj fc Chieffil (I2006D 1. 

For this study, we followed basically the same procedure 
as for the ^®A1 one. We analyzed the same dataset in 
the 1170-1176 and 1330-1336 keV bands, with the same 
templates to estimate the ®°Fe flux. The background de¬ 
termination method lSection l2.4.1l) leads to the identifi¬ 
cation of 2900 and 2961 time-segments respectively. Note 
that the contribution of the ^®A1 line through its inter¬ 
actions with the detectors (Compton effect, diffusion,...) 
has to be taken into account for the ®°Fe analysis. Being 
a strong line, ^®A1 photons can interact inside the cam¬ 
era through the Compton effect and diffusion (the instru¬ 
ment response is non-diagonal) and produce a continuum 
from 20 keV to 1.8 MeV in the data-space. This emis¬ 
sion component is thus convolved with the instrument 
response to predict the corresponding counts in the data- 
space, and the predicted counts in the 1170-1176 and the 
1330-1336 keV bands are included in the sky model dur¬ 
ing the ®°Fe fluxes analysis. Moreover, the contribution 
from the diffuse continuum has been taken into account 
assum ing the power law m odel determined in a previous 
work (iBouchet et al.l[2Mlll . Note that these effects are 
very small, well below the statistical error bars. 

For all the tested distributions, the ®°Fe isotope lines 
are detected at a level of 2cr and 3a respectively, in the 
1170-1176 and 1330-1336 keV bands (Table (H). Their 
global mean-flux in the inner Galaxy is about 4x 10“® 
ph cm“^ s“^. Systematics in the flux determination of 
these large-scale structures due to the background pat¬ 
tern determination method is below ^25% for the ^®A1 
and around ^30% for the ®°Fe lines (fluxes obtained with 
the fixed-pattern are systematically higher). 

As a conclusion, even though poorly constrained, the 
®°Fe to A1 ratio obtained during this analysis is around 
0.14 , which agre es with valu e s prev iously obtained by 
iSmithl (j2004t l and iWang et al.l (120071 1. 

4. DISCUSSION AND CONCLUSION 

For more than a decade, the reference map for the 
spatial morphology of the A1 emission was provided by 
COMPTEL. Several maps have been published, along 
the CGRO mission, with more and more data, but also 
different analysis methods. Consequently, while the ex¬ 
tended morphology of the emission is assessed and results 
compatible, some local features not always appearing, de- 
pending on t he m e thod. Indeed, the firs t MEM images 
IjDiehl et al.l 119951 : lOberlack et al.l Il997tl exhibit many 
low-intensity structures. A small number of them remain 
in the latest MEM image (iPliischke et al.l [200111 . giving 
some confidence in their reliability. In parallel, the chief 
features of these MEM images were confir med with a 
MREM algorithm (|Knddlseder et aLlll999bfl . This algo¬ 
rithm is based on an iterative expectation maximization 
scheme and a wavelet filtering algorithm. This wavelet 
filter suppresses the low-significance features, which are 
potential artifacts, by applying a user-adjustable thresh¬ 
old. This aims to pr oduce the smoothest ima ge con¬ 
sistent with the data (jKnddlseder et al.l Il999b[l . Note 
that the early MREM map built from SPI data by 
iKnodlseder et al.l (j2006l l was too rough to bring any ad¬ 
ditional information compared to COMPTEL ones. 

Now that more than 10 years of SPI data are avail¬ 
able, our main goal was to refine the COMPTEL view 
of the ^®A1 line emission at 1.8 MeV and to investi¬ 


gate the related ®°Fe lines around 1.2 and 1.3 MeV. To 
study their spatial morphology, we have developed spe¬ 
cific tools. We have performed two kinds of image recon¬ 
structions. The first one is based on existing maps, at 
various wavelengths, which are used as a sky model, con¬ 
volved with the instrument response and compared to the 
data (model-fitting). The second one consists in direct 
sky reconstruction from the data, implying an inverse- 
problem. While the first method reveals the global large 
scale morphology of the emission, the latter allows us to 
look for small scale structures like local regions of ^®A1 
production. 

The ^®A1 line is detected at ~ 20cr, to compare to the 
^ 30cr obtained wi th the analysis of COMPTEL data 
obtained in 5 years (|Kn6dlseder et al.lll999all . 

In addition, comparing t he MREM COMPTEL image 
(jKnddlseder et al.l Il999b[ l and SPI one, we report that 
a number of structures appear in both analyzes. Note 
that the SPI sky exposition is non-uniform, mostly con¬ 
centrated along the Galactic plane and differs from the 
COMPTEL one. Thus, SPI is more sensitive in the 
Galactic plane than at high-latitudes, with a point-source 
(3 a) sensitivity of 1.4x 10“® ph cm“^ s“^ in the Galactic 
Center region to be co mpared to the (3 . 8 to 1 .4x 10ph 
cm“^ s“^ reported bv IPliischke et all (|2001h for COMP¬ 
TEL. 

4.1. Data analysis : issues and solutions 

A first point to mention is that, for SPI as for 
COMPTEL, the background treatment is a tricky 
issue. In SPI data, to disentangle the signal and the 
instrumental background contributions, we rely on the 
ability of the instrument to measure simultaneously 
both of them, due to the properties of the coded-mask 
aperture imaging system. Moreover, the evolution of 
the background intensity with time can be determined 
with a segmentation code developed specifically to take 
into account the background variation. The major 
advantage is that it strongly reduces the number of 
parameters related to the background. In a second 
step, we determine the background detector pattern by 
assuming that high-latitude exposures constitute a good 
approximation of “empty-fields”. However, they may 
contain signal from the diffuse emission (Section 12.4.21) 
We have thus implemented the possibility to determine 
the background pattern during the data-reduction 
process, assuming that the signal is sufficiently well 
known. The two approaches for background pattern 
determination have been compared to assess the robust¬ 
ness of the results, leading us to essentially the same 
conclusions regarding the ^®A1 emission characteristics. 

4.2. ^^Al and large-scale emission 

In the model fitting analysis based on other wave¬ 
length maps, our data do not allow us to distinguish a 
unique preferred template since several ones lead to sim- 
ilar likelihood parameter values. From COMPTEL data, 
iKnddlseder et al.l (|1999al) had concluded that the best 
tracers were DIRBE 240/r and 53 GHz free-free maps. We 
agree that the FIR maps with wavelengths 60 to 240 /x 
or 53 GHz (free-free and dust) appear statistically as the 
best estimates of the ^®A1 emission global morphology 
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but CO, NIR A[4.9/x] and A[3.5/i] extinction-corrected 
maps have also to be considered with only a slightly lesser 
degree of confidence. In addition, “low-contrast” maps 
observed at 25p, and 12^ provide an equally good de¬ 
scription of the emission. 

Indeed, our results confirm that the Al emission follows 
more or less the distribution of the ext reme Population / , 
the most massive stars in the Galaxy (iDiehl et al.lll99^ . 
It is known that the massive stars, supernovae and no¬ 
vae produce the long-lived isotopes ^®A1 and ®°Fe with 
half-lives of 0.7 and 2.6 My. On the another hand, large 
amounts of dust/grains condense in core collapse SNe 
ejecta while massive star-supernovae are major dust fac¬ 
tories, therefore dust FIR and MIR maps are expected 
to correlate with the corresponding emissions. 
Quantitatively, the flux extracted through the model fit¬ 
ting method for the inner Galaxy {\l\ < 30°,|6| < 10°) is 
found to be around 3.3x 10“'^ph cm“^ s“^ while it does 
not depend much on the sky model (particularly if we re¬ 
strict ourselves to the preferred ones). This flux is consis¬ 
tent with earlier measurements of both INTEGRAL SPI 
and COMPTEL instruments and corresponds to a to¬ 
tal ^®A1 mass contai n ed in our Galaxy of ^3M0,(bu t, 
see also IDiehl et all (I201C) and iMartin et al.l (1200911 '). 
However. iChurazov et al.l ( 20101) using a ’’light-bucket” 
method obtain a flux of 4 x 10“^ ph cm“^ s“^ for the 
region delimited by |/| < 40°,|6| < 40°. For the same re¬ 
gion, we obtain a flux of ^ 3.9 x lO^'^ph cm“^ s“^ with 
the most probable spatial morphologies. 

Another observable quantity, related to ^®A1 production, 
is the ratio between ®°Fe and Al line fluxes in the inner 
Galaxy. From a theoretical point of view, the ®°Fe to Al 
flux ratio provides a test for stellar models, as predictions 
of the yields of massive stars depend strongly on the 
prescripti on of nuclear rates, stel l ar winds, mixing and 
rotat ion (['Wooslev fc HegeJ l2007t iTur. Heeer fc AustinI 
1201011 . With our analysis, the ®°Fe mean flux in the inner 
Galaxy is about 4x 10“® ph cm“^ s“^, leading to a ®°Fe 
to ^®A1 ratio between 0.12 and 0.15. Considering the 
uncertainties still affect ing the models, thi s results con¬ 
firms those reported by iHarris et al.l (|2005ll ; iWang et al.l 
(|2007ll . and can be used to reject some hypotheses, but 
not yet to definitively discriminate the good ones. 

4.3. Imaging results 

To go further in the details of the emission distribu¬ 
tion, the flux extraction must be independent of any 
template, and rely on a direct imaging reconstruction of 
the ^®A1 emission. To realize it, we chose the MEM tool 
because of its ability to pro cess high-resolution ima ges. 
Our MEM code is based on iSkilling fc Brvaril (|1984[ 1 al¬ 
gorithm. The main improvement we implemented to the 
basic algorithm is the possibility to refine the background 
determination during the iterative computation of the so¬ 
lution. 

The resulting SPI image presented in Figure [5] resem¬ 
bles the MEM COMPTEL images in term of angular 
resolution and details. With a resolution fixed to ^6°, 
it gives essentially the same information as the previ¬ 
ous method, i.e. the emission is mainly confined in the 
inner-Galaxy disk (with |6| < 7°). But, it also suggests 
the presence of extended (a few degrees) emitting ar¬ 
eas and makes it quite instructive to compare the main 


features with those observed in COMPTEL maps: Sev¬ 
eral excess es appear in bo t h SP I and COMPTEL data 
analyzes of lPliischke et al.1 (120011) . which reinforces their 
reliability see (Table [2|). 

In the Cygnus region, Cyg OB2 cluster detection has 
been reported in COMPTEL and SPI data (jMartin et al.l 
I2009D . We find at the same position a flux of 4.1 — 4.5 ± 


1.5x 10 


-5 


ph cm 
1-5 


-2 S -1 


compatible with the value of 

-2 „-i - 


3.9 ± l.lx 10 ° ph cm ^ s ^ reported bv IMartin et al.l 
(I2009D and 3.7 ± 1.1) x 10“^ ph cm“^ s“^ 


^ reported by 

iPliischke et al.l (|2000l) . However, a more complete anal¬ 
ysis based on the excesses visible in the image reveals 
several spots in the Cygnus region and around (Table 2), 
suggesting a complex structure of this area. 

The Vela and Carina regions are detected, in our work, 
at a low significance level (2cr). For Vela, t he measured 
flux is comparable to the value reported by IDiehl et al.l 
(|1995D . For Carina, th e flux is comparable to the 
estimation of IVoss et al.l (I2012D using SPI data and 
iKnodlseder et al.l ~ 1996D using COMPTEL data. 

In the Sco-Cen region, emission, not detected in the 
COMPTEL data, was reported in the SPI data at a 
high signifi c ance level (6 ± lx 10“® ph cm“^ s“^) by 
IDiehl et al.l (|2010D with a somewhat different analysis. 
Using their location and spatial extension, we find a flux 
of 4.2 ± 1.6x 10“^ ph cm“^ s“^ in the most favorable 
case (the flat H/ map used to model the diffuse large- 
scale emission). In our image, the excess appears shifted 
and less extended (Table 2). Moreover, the measured 
flux depends strongly on the template used to model the 
large-scale emission of the ^®A1 fine. This lead us to 
consider that this result requires additional work to be 
confirmed. 

On the other hand, we pick up significant and robust ex¬ 
cesses, not explicitly reported by the COMPTEL team. 
Two of them, rather extended, are seen in the Tau¬ 
rus/Anticenter region at (l,b)~ (161°,-3°) and (l,b)~ 
(149°, 8°) with fluxes of 4 — 6 ± 2 and 8 — 9 ± 3x 
10“® ph cm“^ s“^. It is not excluded that they be¬ 
long to a same structure, leading to a ^4tT detection. 
Furthermore, they can be correlated with a weak feature 
aro und (1, b) -(160°,0°) in t he COMPTEL MREM im¬ 
age (IKnodlseder et al.lll999bD . A third feature appears at 
high-latitude ((7 ± 2)x 10“^ ph cm“^ s“^). Obviously, 
we can not exclude that this spot is due to statistical 
data noise, but its stability against the various analysis 
methods supports its reliability. 

We have to point out that the present analysis has been 
optimized for large scale emissions (handling the all-sky 
data set) while any ’’local” analysis requires a dedicated 
procedure. Briefly, in this case, we have to use the expo¬ 
sures whose pointing direction are not too far from the 
source or region of interest (less than 12°) and to refine 
the sky mo d el, to optimize th e signal-to-noise ra t io. 
iChen et ^ (|1996D (see also iPrantzos fc Diehll (119951 )1 
have proposed an interpretation of the COMPTEL ^®A1 
emission by linking enhanced emission spots to the 
Galactic spiral arms. This relation has been con¬ 
firmed with the observation of the line shift along 
the Galactic plane in the SPI data (|Diehl et al.l 120061 : 
iKretschmer et al.ll2013D . Some of the identified regions 
match excesses we point out in our analysis, reinforcing 
the reliability of the emissions and supporting the pro- 
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posed scenario (labeled with ’ * ’ in Table 2). 

To conclude, the SPI instrument gives us one of the 
rare opportunities to get pieces of information about the 
spatial distribution of ^®A1 line and the nucleosynthe¬ 
sis process. Fifteen years after the first results obtained 
by COMPTEL, we confirm the chief features as well as 
some of the low-intensity structures reported by this in¬ 
strument through direct-imaging of the sky and tem¬ 
plates comparison. While we cannot hope to increase 
significantly the amount of data recorded by INTEGRAL 
(more than 10 years of observation have been included in 
the presented analysis), some improvements in the data 
analysis appear achievable to still refine our results. In 
particular, an accurate modeling of the response outside 
the field of view is a prime objective for studying emis¬ 
sions above ~ 1 MeV together with the analysis of the 
multiple events, which contain ~20% of the y-ray pho¬ 
tons. 
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APPENDIX 

THE MAXIMUM ENTROPY ALGORITHM 


The Maximum-Entropy algorithm aims to maximize the following function Q (Equation 4): 


Qif) = 


C{f) 

2 


(Al) 


where a is a regularization parameter. Eor an image containing N pixels, the following entropy function H is used 


f- 

H{f,m) = '^{f^ -rrii- f^log—) 
^ rui 


(A2) 


where rrii is a model image assigned to pixel which expresses our prior knowledge about the sky intensities fi.The 
sky intensities and their model are positive quantities. C{ f) measures the discrepancy between the measured data d 
and the reconstructed model of the data. A single statistical constraint is generally used and for data with a Gaussian 
noise, C is the chi-square function: 


C = 


M 

E 


(E^i R^jfj -b^- 


(A3) 


In this expression M is the number of measured data-points rij, represents the noise level and bi the instrumental 
background for each data-point. i? is a matrix of size M x N, which represents the response of the instrument. To 
simplify the presentation of the algorithm, the values bi a re su pposed to be known and provided by the user. 

The resulting map / obtained by maximizing Equation lAll will be unique since the surfaces H = constant and 
C = constant are both convex. For every fixed constraint level, Caim, 

C{f) = Ca^m (A4) 

defines an ellipsoidal hyper-surface of radius Caim in t he N -dimensional image-space. For this surface there is only one 
tangent point / with a certain (maximal for Equation I A4() entr opy hyper-surface for which the gradient VQ vanishes. 
Therefore, f is the maximum entropy solution for Equation IA4I given by the formula: 

VQ = aVH - ivC = 0 (A5) 


Convergence method 

An iterative procedure is needed to solve the set of non-linear equations lEouation IA5I) . It can be the following: 
Starting from the point of absolute maximum-entropy solution { = m, a = oo ), one finds at each iteration, 

the maximum entropy point / , which fulfills the condition of Equation IA5I for certain ellipsoidal hyper-surfaces 
(Ecuation IA4I) with the condition: 

Cif + 6f)<Cif) (A 6 ) 

The new solution is updated by incrementing the previous solution with a vector Sf. The solution is forced to follow 
the maximum-entropy trajectory by using the appropriate value of a defined by Equation IA5I 

As the solution is updated at each iteration, the approximation of Q, through its quadratic expansion, should remain 
accurate. Maximizing Q(/) subject to |<5/p < Iq for some small value of fo, between successive iterations, fulfills the 
requirement. In addition, the tendency of the search-direction algorithm to produce negative valu es of / is limited. 
Howe ver such distance limit tends to slow the attainment of high-values of/. To overcome this defect, (|Skilling &: BrvanI 
1198411 (hereafter SB84) suggest to use, instead of the squared length of the increment: 

I" = < ^0 (A7) 

i 

Using a distance in this form is equivalent to put a metric with —WH onto the image-space. 

= 1 /fii = j) and = 0(f ^ j) (A 8 ) 

The upper index denotes a contra-variant vector. 

However, the problem is a N-dimensional problem and can be difficult to solve for a large number of unknowns N with 
a Newton-Raphson procedure. 

It is one of the reasons why we use the algorithm proposed by SB84. It consists in choosing, at each iteration, three 
search-directions ei, 62 and 63 derived from the gradients VC and Vi/, as well as from the matrix of curvature of C, 
VVC. In this way, the N-dimensional problem is reduced to a three-dimensional problem. The three search-directions, 
by using the metric given by Equation IA 8 [ are 

ei = /(Vi/), 62 = /(VC), 63 = |Vi/|-V(VVC)/(Vi/) - |VC|-V(VVC)/(VC) (A9) 
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The entropy metric is used to define the lengths 




^df 

The solution is updated by adding the increment vector as follows 

3 


/ = / + 5/ = /+E 




(AlO) 


)i=i 


The values to be found are the lengths ( ^=1,2,3) along the search directions e^. For this purpose, H and C are 
modeled by their Taylor 2 -order expansions H and C at the current solution / 


H{x) = Hq + H^x^ - , C{x) = Cq + C^x>^ - l/2M^^x^x''] 


where 


= VH, • VC, 


: 9t^l^ — • WC • 


After simultaneous diagonalisation of the quadratic models of H and C, in a common basis, formed by the eigenvectors 
of the matrix computed on the basis formed by the eigenvectors of the matrix one obtains: 

IL{x') — Hq -\- Hfj^x^ l/2x^x^, C(x) — Cq -t- C^x^ l/2y^^jx^x^, I — x^x^ 

where 7 (^) are the eigenvalues of This is a low-dimensional problem and standard algorithms can perform these 

operations. According to the maximum-entropy trajectory (VQ = 0 as defined in Equation IA5I) . the step-lengths are 


{aS^ - C^) 


" {7(m)+«) 

SB84 introduce explicitly the distance constraint into the optimization process.The quadratic model of Q becomes 


(All) 


Then, the step-lengths x^ are given by 


Q = ail — (7/2 -|- pP/2 with p > 0 


- Cf, 

P + 7(/x) + “ 


(A12) 


(A13) 


In SB84, p can be interpreted as an artificial increase of each eigenvalue 7 (^) of C. The equivalent quadratic model of 
C is 

Cp{x) = Co + CpX^ + 1/2(7(^) -b p)xpXp (A14) 

Control of the algorithm 

The control of the convergence of the algorithm is performed directly on the constraint C, and not on the Lagrange 
multiplier a. The aim is to maximize S over C = Caim- The minimum reachable value of C is 

Crmn = ^0 - - V -^ 

SB84 recommend to have a more modest aim Caim, for example 

Caim — IXiaxi^/i^Cmin “b l/3Go, Caim') 

C{x) is an increasing function of a. The simultaneous control of C{x) and is done through the value of p, C{x) is 
an increasing function and P a decreasing function of p. By adjusting the values of a and p, one can reach the aimed 
result. The procedure is detailed in SB84. 


Stopping the algorithm 


In the ” historical” version of the maximum-entropy, this process is repeated until C reaches the stopping value M. 
SB84 suggest to reach at least a value lower than the largest acceptable value at 99 % confidence, about M + 3.39-\/M 
where M is the number of observations. 

Equation IA5I forces the gradients VH and VC to be parallel. Then, the parameter a can be interpreted as the ratio 
of their lengths. In addition, the algorithm is always checked by measuring the degree of non-parallelism between VH 
and VC through the value of 


TEST = 


1 VH 
r \VH\ 



(A15) 
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The value is zero for a true maximum entropy image. SB84 indicate that reaching TEST < 0.1 or so at the correct 
value of C, demonstrates that the unique maximum-entropy reconstruction has been attained. 

In our application, we always reached C < M and TEST <0.1. The function C always decreases when the iteration 
number increases. From a number of iterations, the value of TEST starts to decrease in a monotonic way and reaches a 
minimum value. From this point, its value becomes almost stable and in the worst case increases slightly whereas there 
is essentially no progress in the decrease of the value of C. In addition, we compare the function C to its quadratic 
approximation to ensure its validity. To achieve and maintain this constraint during the iterations, the length of the 
increment ('Eauation lA7l) is kept small enough, this results in a larger iteration number. In general, 30 to 50 iterations 
are required (less if the increment constraint is higher or relaxed at any given number of iterations, but at the expense 
of an inaccurate 2nd-order approximation of C). 


Our application 

Different choices of C (f) are possible (jBrvan fc Skilling1ll980H We u s ed a modified version of the statistic, 
accurate for low numbers of counts, following the prescription of iMight^ (|1999tl . The image model is computed at 
each iteration as 

rrii = cste = exp [—-) 

Si=l fi 

SIMULATIONS OF IMAGING RECONSTRUCTION 

To simulate synthetic data, we rely on the model-fitting analysis. A given template map is used to model the 
distribution of the ^®A1 line over the Galaxy (Section 3.2). Expected counts are obtained by convolving this sky model 
with the instrument response and then adding the background model. The intensity of the map and the intensity 
variations of the background are the parameters which are adjusted to the recorded data. 

Our simulations start from the predicted counts. Poisson statistical fluctuations are added to them to build the 
simulated data, which are in turn analyzed similarly as the real data. We simulated data from the MREM template, 
which represents a smooth input map, and EGRET template, which is a more structured. The input maps and their 
reconstruction are displayed on Figure |51 We present two image reconstructions which differ in the resolution of the 
final image, forced to the values of 3° and 6° (FWFIM). In all cases, the fluxes simulated and measured in the inner 
Galaxy are recovered within 5% for both MREM and EGRET simulated maps. From these simulations, we conclude 
that statistical noise produces structures whose intensity can reach about one-tenth of the image maximum intensity. 
Above this value, there are no particular structures created during the image reconstruction process at any particular 
position. The longitude profiles of the reconstructed images are compared to the simulated images on Figure [51 
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Figure 1. Mean count rates recorded on the camera per detector versus time in the 1805-1813 keV energy band. In gray, the count 
rates and associated error bars (1 point per exposure). In black, the count rates averaged per revolution. The total number of exposures 
is 76 789. The hole around days 2500 corresponds to exposures with no PSD information (Section [JJ. 



Figure 2. A sky model, consisting of a synthetic map (here the 60/i map. Section l3.2fl . plus a few sources, is fitted to the 1805-1813 
keV data (fixed-pattern method). The evolution of the reduced {x%l^of) (solid-line, open black-circles) and the x% ~ (dashed-line, 
open red-diamonds) are shown as a function of the assumed background timescale (in hours). In insert are displayed the x?, (dotted-green 
line) and the number of parameters to be determined (solid black line). The filled black-circle, red-diamond and (in insert) green triangle 
and black-square show respectively the the x?^ ~ ^^e number of parameters obtained with the background 

segmentation method. The infinity symbol (oo) on the x-axis corresponds to a constant background for the whole dataset. 
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Figure 3. Background intensity evolution in revolution 392. The background is modeled with 3 time-segments (dotted purple-line) to 
give a x^/do/ of 1.02 or with 5 time-segments (red line) for a jdoj of 0.87. 
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Figure 4. Flux of the sky obtained with the fixed-pattern (open red-circles) and fitted-pattern (filled black-circles) background method 
as a function of the assumed background variation timescale (here using a 60^ map). The values obtained with the optimized number of 
background parameters are shown with the squares (open for fixed and filled for fitted-patter). 





































































Galactic ^®A1 emission with INTEGRAL SPI 


15 



1.03x10" 1.33x10'® 2.56x10'® 

ph.cm“®.s“\sr"^ 


Figure 5. Image of the line (1805-1813 keV). The ima ge is b uilt with a resolution (ICF FWHM) of 6°. For this image reconstruction, 
the background pattern is adjusted as explained in Section 12.5.31 The contours are extracted from the 3° resolution image. In units of 
xl0~^ ph cm“^ s“^ sr“^, they correspond to 0.54, 1.1, and 2.7. Identified regions, from left to right: Perseus region (105° < 1 < 170°) 
(Taurus clouds), Cygnus/Cepheus region (75° < 1 < 100°), the inner Galaxy (-30° < 1 <30°, -10° < b <10°), Carina (1=286°, b=l°) and 
the Vela region (260°< 1 < 270°). At mid-latitude, the Sco/Cen region (300° < 1 < 360°, 8 °< b < 30°). 



map 


Figure 6. The relative chi-square variation ^) versus assumed template to model the distribution of the ^®A1 line. ^) the 

from which the value of the best fitted template is subtracted. The best template is the 100/i template for the fitted-pattern method 
(;^2 _]^258778.1 for 1261797 dof) and the 25// template for the fixed-pattern method (x|^=1266968.4 for 1262208 dof). Terms: sync., dust 
and free-free, are abbreviations for 53 GHz synchrotron, dust and free-free maps described in Table ^ MEM and MREM indicate the 
COMPTEL maps, A[3.5/i] and A[4.9/i]corrected NIR extinction map. The red curve is for the fitted-pattern and green for fixed-pattern 
method. The maps are ordered following their contrast defined as the ratio of the flux enclosed in the region |/| < 150°, \b\ < 15° to the 
total flux. 
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map 


Figure 7. Flux in the inner Galaxy as a function of the map used to model the distribution of the ^®A1 line. The labels are the same 
as in Figure [6| The dashed black-curve is the total flux in the Galaxy (fitted-pattern), scaled by a factor 0.2. The dotted lines (red for 
fitted-pattern and green for fixed-pattern) are the fluxes obtained, if an isotropic emission (possibly extra-galactic) estimated by using the 
map emission at |6| >40°, is subtracted from each template. The labels are the same as in Figure[6] 


Table 1 

Maps used as templates. 



Tracer 

Mechanism 

fMIR 12 and 25/i 

warm dust (T'^250 and ^^120 K) 
/AGBs 

dust nano-grains and PAHS 
heated to high temperature 

t53 GHz sync 
“ COMPTEL-MEM 

cosmic-rays/magnetic field 

Synchroton 

IC 

Inverse-Compton from GeV 
electrons on the CMB and ISRF 

Inverse-Compton 

fNIR 1.25, 4.9 /.t 

stars (K and M giants) 

star light 

(21 cm) 

H hyperfine transition 

Neutral hydrogen 

fNIR 2.2/7 

stars (K and M giants) 

star light 

fEGRET 

Interstellar gas/cosmic-rays 

nuclear interactions 

fNIR 3.5 fj, 

stars (K and M giants) 

star light 

t53 GHz free-free 

ionized gas 

Free-free 

t53 GHz dust 

dust 

Thermal dust 

fFIR 100, 140 and 240/i 

“COMPTEL-MREM 

warm dust (T ~30, 

-^21 and --12 K) 

microns sized dust emitting in thermal 
equilibrium with the heating ISRF 

fFIR 60/t 

warm dust (T'^hO K) 

microns sized dust emitting in thermal 
equilibrium with the heating ISRF 

'^CO 

GO rotational transition 

Molecular gas / young stars 

‘^NIR extinction-corrected map 

3.5, 4.9 (hereafter A[3.5/i and A[4.9/i] ) 

stars (K and M giants) 

star light 


Note. — The maps are ordered in ascending “contrast”, defined here as the ratio between the fraction of the emission 
enclosed in the region |/| < 150°, \b\ < 15° to that of the whole sky. The value of the ratio varies from 0.4 (25^) to nearly 
1 (A[4.9 m]). 

t available at http://lambda.gsfc.nasa.gov or http://heasarc.gsfc.nasa.gov. 

The NIR, MIR and I' lK maps are the UURE/DIRBE Zodi-Subtracted Mission Average (ZSMA) maps, from which zodi¬ 
acal light contribution has been subtracted. Our 53 Ghz synchrotron presents some inaccuracies and is used for indicative 
purpose. “ The maximum-entropy (MEM) and Multiresolution Regularized Expecta tion Maximization (M REM) all-sky 
image of the Galactic 1809 keV line emission observed with COMPTEL over 9 years flPliischke et al.ll2001L and references 
therein). teickev fc LockmanI II1990I) . ^Dame, Hartmann &: Thaddeuj 11200 U ) with central peak removed fSection 12.311 . 
^The NIR extinction maps fSection l2.3ll . 


















Galactic emission with INTEGRAL SPI 


17 



o 

■o 

3 


iS 

o 

'i3 

O 

jS 

<0 


C5 


EGRET-6 deg. 





(a) Image template (shown at 6° FWHM) 



(b) Reconstruction 6° FWHM 



Figure 8. Simulated and reconstructed maps. The image contains 49152 pixels (pixel size is 0.9°) ordered following the Healpix scheme. 
(Top) Simulated MREM (initially 3.8°) and EGRET (initially 2°) displayed using an angular resolution of 6° FWHM. (Middle) The 
reconstructed MREM and EGRET fixing the output map resolution to 6° FWHM and (Bottom) to 3° FWHM. Images are displayed using 
the DS9 software (http://ds9.si.edu/site/Download.html). The images are scaled between their maximum intensity fmax and fmaxT^^- 
The scale is in square root unit ot the tlux (ds9 Scale Square Root menu). (Top) contours correspond to /maa:/2, /maj;/4, fmaxA^ and 
/ma3:/l6* (Middle) fmax fmaxj^ and /ma3:/l0. (Bottom) fmaxj^ and fmaxN‘^- 
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(a) EGRET simulated template 



(b) MREM simulated template 


Figure 9. Longitude profiles of the reconstructed images compared to the simulated images (See Figure [8] 


Table 2 

^®A1 emission sites and flux in units x 10“^ ph cm~^ s“^. 


Source 

name 

Position 

1°, b° 

Spatial 

morphology 

^®A1 large-scale morphology 
IC map A[4.9fj] 



Cygnus region 



Cyg OB2 cluster‘d 

81,-1 

Gaussian, (t=3 

4.5 ± 1.5 

4.1 ± 1.5 


Mutiple 

sources - most significant excesses 



t81,-l 

Gaussian, cr=3 

4.3 ± 1.6 

3.8 ± 1.6 


100, 6 

Point-source 

2.9 ± 1.1 

2.7 ± 1.1 


+ 64, 1 

Point-source 

2.2 ± 1.1 

2.0 ± 1.1 

Sco-Cen“’+ 

350, 20 

disk r=10° 

0.4 ± 1.4 

1.9 ± 1.4 


360, 16 

disk r=5° 

0.5 ± 0.9 

1.1 ± 0.9 



Orion/Eridanus region 



Orion/Eridanus^ 

150 < 1 < 210, -30 < b < 5 

1.2 ± 3.2 

0.5 ± 3.1 

Carina^ 

*287,0 

disk r=3° 

3.1 ± 1.5 

2.8 ± 1.5 

Vela 

*267,-1 

disk r=4° 

3.6 ± 1.8 

3.3 ± 1.8 

Taurus^ 

*161,-3 

Disk r=5° 

5.6 ± 2.1 

4.4 ± 2.1 


149, 8 

Disk r=ll° 

8.5 ± 2.9 

7.9 ± 2.9 



Other significant excesses 




226, 76 

Gaussian, (T=3° 

7.2 ± 1.8 

6.8 ± 1.8 


Known high-energy emitting sources 



Crab 

185,-6 

Point-source 

1.2 ± 1.3 

1.1 ± 1.3 

Cyg X-1 

71, 3 

Point-source 

0.2 ± 1.1 

0.2 ± 1.1 

PSR B1509-58 

*320,-1.2 

Point-source 

2.8 ± 0.9 

2.8 ± 0.9 

Note. — The first 

column contains the name of the source if it 

is known. The 

second column 


the position of the source in galactic coordinates and the third column the source spatial morphology; 
a point (point-source), an axi-symmetric Gaussian {cr indicated) or a disk (radius indicated)). The 
position and spatial morphology in bold correspond to known sources or are based on published works. 
Next columns give the source fluxes obtained for two template maps used to model the large-scale 
distribution of ^®A1 over the Galaxy; IC (low-contrast map) and A[4.9/i] (high-contrast map). A 
model, comprising the ^®A1 line distribution, the sources and the background, is adjusted to the data. 
The uncertainties are obtained from the model-fitting analysis by using the curvature matrix of the 
likelihood function and terms associated with the variance of the solu tion. The differen t lab els identify 
excesses, which have been already mentioned in the literature by: “ iPiehl et ahlf^lOll - ^ lIVoss et al.l 
^ iPiehl et al.lfTo^l - ^ JMartin et al.l|200^ - ^ JPiehl et al.ll2002ll . 

* Possible association with the spiral arms iChen et al.llT996D - ^ Suspected or visible in COMPTEL 
^®A1 image iOberlack et aOl 19971: iPliischke et al.ll200in . '*'See Section l3.3l 


Table 3 

Radioactive line fluxes in the inner Galaxy in units of xl0~'^ ph cm~^ s~^ . 



25/1 

COMPTEL-A\ 

2i0/i 

COMPTEL-MKEM 

60/1 

CO 

t A[4.9 m] 

26 Al 

3.68 ±0.21 

3.41 ±0.20 

3.15 ± 0.18 

3.24 ±0.19 

3.11 ±0.18 

3.05 ±0.18 

2.92 ± 0.17 

60 Fe 

0.47 ±0.16 

0.45 ±0.15 

0.38 ±0.14 

0.38 ±0.14 

0.38 ±0.13 

0.41 ±0.13 

0.42 ± 0.13 

60Fe/26Al 

0.13 ±0.05 

0.13 ±0.05 

0.12 ±0.05 

0.12 ±0.05 

0.12 ±0.05 

0.14 ±0.05 

0.15 ±0.06 


Note. — Fluxes are given for the inner galaxy {\l\ < 30° and |b| < 10°). f "the map is corrected from reddening and fluxes are 
essentially zero for |6| > 20°. ®*^Fe is the average flux of the 1173 and 1333 keV lines (1170-1176 keV and 1330-1336 keV bands) 
and ^®A1 the flux of the 1809 keV line in the 1805-1813 keV band. The fluxes are obtained using the “fitted-pattern” method. 
































